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Abstract 

In this paper we study three self-similar, long-range dependence, Gaussian processes. The 
first one, with covariance 

rsAt 



u^lit-u)" + {s-u)'']du, 







parameters a > —1, — 1 < 6 < 1, |6| < 1 + a, corresponds to fractional Brownian motion for 
a = 0,— 1<6<1. The second one, with covariance 



i2-h)(^s'^ + t'^-l[{s + tf + \s~t\'^]y 



parameter < ft. < 4, corresponds to sub-fractional Brownian motion for < h < 2. The 
third one, with covariance 

log s + t^ log t^^[{s + t)^ log(s + t) + {s - \og\s - t\ 

is related to the second one. These processes come from occupation time fluctuations of 
certain particle systems for some values of the parameters. 
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1 Introduction 



Some types of self-similar, long-range dependence, centered Gaussian processes have been ob- 
tained from occupation time fluctuation limits of certain particle systems (they represent the 
temporal structures of the limits). These processes have parameters whose ranges of possible 
values are determined by properties of the particle systems. Since they are interesting processes 
in themselves, it is natural to ask if the ranges of values of the parameters can be extended 
and what the maximal ranges are, what properties the processes possess, and if the extended 
processes may also be related to some stochastic models. In this paper we will consider the 
first two questions for some extensions of fractional Brownian motion (fBm) and sub-fractional 
Brownian motion (sfBm). The third question, which is of a different kind (weak convergence of 
processes related to the particle systems), is discussed in [6|. 

In the remainder of the Introduction we will discuss the motivations behind the present 
article in connection with particle systems. In the following sections we will study the processes; 
most of the contents in those sections can be read independently of this discussion. 

We will study a centered Gaussian process {Ct)t>o with covariance Q{s,t) = E^sCt, s,t > 0, 
given by 

Q{s,t)= u^iit - u)^ + (s - u)^]du. (1.1) 
Jo 

For a = 0, (jl.ip becomes 

+ (1.2) 

which for — 1 < 6 < 1 corresponds to the covariance of the well-known fBm with Hurst parameter 
^(6+1), Brownian motion (Bm) for 6 = 0. We call the process ^ with covariance (II. ip a weighted 
fractional Brownian motion (wfBm), refering to the weight function u*^ in the integral. Thus, 
wfBm's are a family of processes which extend fBm's. 

The process ^ appeared in [5] in a limit of occupation time fluctuations of a system of 
independent particles moving in M*^ according a symmetric a-stable Levy process, < a < 2, 
started from an inhomogeneous Poisson configuration with intensity measure dx/{l + |x|'''), 
0<7<d = l<a, a = — 7/a, 6 = 1 — 1/a, the ranges of values of o and 6 being —1 < 
a < and < 6 < 1 + a. The homogeneous case (7 = 0, fBm) was studied in [21 [3]. The 
process ^ also appears in [6j in a high-density limit of occupation time fluctuations of the above 
mentioned particle system, where the initial Poisson configuration has finite intensity measure, 
with d = 1 < a, a = —1/a, 6 = 1 — 1/a. From these results it follows that for —1 < a < 
and any < 6 < (1 -|- a) A (1/2), the kernel (II. Ih is positive-definite, because it is a limit of 
covariances. The question that interests us here is what are the maximal ranges of values of 
the parameters a and 6 for which (jl.ip is a positive-definite kernel (the conditions a > —1 and 
6 > — 1 are obviously necessary). We will show that the maximal ranges are a>— 1,— 1<6<1, 
and |6| < 1 + a, and we will study some properties of the Gaussian process with covariance 
(jl.ip . These properties, which are analogous to those of fBm, are self-similarity, path continuity, 
behavior of increments, long-range dependence, non-semimartingale, and others. Since wfBm 
widens the scope of behaviour of fBm, perhaps it may be useful in some applications. 

A different extension of fBm is the hi- fractional Brownian motion (bfBm) introduced in [11] , 
with covariance 

^,2H ^^2H^k _^^_^^2Hk^ (1.3) 

where < < 1 and < < 1, corresponding to fBm with Hurst parameter H for k = 1. 
This process was obtained analytically from the covariance of fBm, not from a specific stochastic 
model, whereas the process with covariance (jl.ip appeared first in connection with a particle 
system. The process bfBm is also studied in [141 119j . With H = 1/2, bfBm can be extended for 
1 < A; < 2, and in this way it plays a role in our next process, which we describe presently. 
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Let (Ci)t>o be a centered Gaussian process with covariance K{s,t) = E^sCti s,t >0, given 

by 

K{s, t) = {2- h) (^s^ + - + tf + \s- tl'^]) , (1.4) 

where < /i < 4, which for h = 1 corresponds to Bm. This process was found in [21 E] 
for 1 < h < 2 (where the covariance was written without the factor (2 — /i), which is irrelevant 
there) in a hmit of occupation time fluctuations of a particle system, which is as the one described 
above in the homogeneous case (7 = 0), and in addition the particles undergo critical binary 
branching, i.e., each particle independently, at an exponentially distributed lifetime, disappears 
with probability 1/2 or is replaced by two particles at the same site with probability 1/2. (See 
also [12j in the setting or superprocesses with a = 2.) In this model d is restricted to a < d < 2a, 
and h = 3 — d/a; hence the range of values of h for the particle model is 1 < /i < 2. However, 
existence of the process ^ for all < /i < 2 follows from the equality in distribution 

{Ct)t>o = c{Ct + C-t)t>o, (1.5) 

where ^ is fBm defined for all —00 < t < 00 and c = I/V2 |2j. The process C with < h < 2 was 
called sub-fractional Brownian motion in [2j because it has properties analogous to those of fBm, 
and it is intermediate between Bm and fBm in the sense that its increments on non-overlapping 
intervals are more weakly correlated and their covariance decays faster than for fBm. Properties 
of sfBm have been studied in [2\ and [17^ I18j. The process (jl.Sp with c = 1/2 is called the even 
part of fBm in [8j, where it appears in a different context. 

We will show that K{s,t) is a covariance also for 2 < h < 4 (degenerate for k = 4). Since 
the coefficient (2 — h) in K{s, t) changes sign from positive to negative as h passes from below to 
above 2, we call the process C for 2 < /i < 4 a negative sub-fractional Brownian motion (nsfBm). 

Note that in contrast to this extension of sfBm, in the case of fBm it is not possible to extend 
the covariance for Hurst parameter > 1 (the kernel is not positive definite). 

It turns out that nsfBm for any 2 < /i < 4 is represented as the integral of the extended 
bfBm referred to above for k = h — 2 (with H = 1/2), i.e., with covariance 

(s + t)^-2 - |s-t|'^-2. (1.6) 

This implies that nsfBm is a semimartingale, as opposed to sfBm, which is not. We will also 
discuss other properties of nsfBm. 

For 2 < h < 5/2, the covariance (jl.4p is a limit of covariances of occupation time fluctuations 
of the branching particle system described above in the homogeneous case for d = 1 < a, with 
h = 3 — d/a (therefore it is a covariance), but this covariance does not correspond to the 
occupation time fluctuation limit of the particle system, because the system becomes extinct 
locally for d < a [lOj. For d = a there is a functional ergodic theorem [15]. On the other hand, 
the high-density occupation time fluctuation limits of the branching particle system obtained in 
[6] yield both sfBm for a < d < 2a and nsfBm for d < a. In this setup, the case d = a leads to 
a centered Gaussian process {r]t)t>o with covariance R{s,t) = Eijsrjt, s,t > 0, given by 

R{s, t) = - (^s^ log s + log t - 1 [(s + tf log(s + t) + {s- tf log |s - (1.7) 

[6], and we will also study this process. Although it has a similar interpretation to that of nsfBm 
regarding the particle system, it is not a semimartingale. 

As we have seen, the processes described above arise in the limits of certain particle systems 
for some values of the parameters. It would be interesting to find particle systems that produce 
these processes also for the remaining values of the parameters. 



We remark that there exist more general models that yield, in the limit, stable, self-similar, 
long-range dependence processes, which in the Gaussian case coincide with the processes above 
for some of their parameter values [21 [3] . 

The present paper was written separately from [5l [6] because it can be read independently. 
The results and the methods of proof in those papers refer to convergence of occupation times 
of particle systems. 

Self-similar processes with long-range dependence have attracted much attention recently for 
their applications and their intrinsic mathematical interest (see e.g. [3116]). Therefore it seems 
worthwhile to investigate new processes of this type, and stochastic models where they arise. 



2 Weighted fractional Brownian motion 
2.1 Results 

Let us recall the funtion Q defined by (11. ip : 

Q{s,t)= u^iit - + {s - u)'']du, s,t>0. (2.1) 
Jo 

We start by discussing conditions under which this function (symmetric, continuous) is the 
covariance of a stochastic process. 

Theorem 2.1 The function Q defined by ^2.1\) is positive- definite if and only if a and b satisfy 
the conditions 

a>-l, -1<6<1, |6|<l + a. (2.2) 

The proof of this theorem will be given in the next subsection. 

Definition 2.2 The centered Gaussian process ^ = {£,t)t>o with covariance function (12. ip satis- 
fying conditions (|2.2p will be called weighted fractional Brownian motion (wfBm) with parameters 
a and b. 

Remark 2.3 (a) It is immediately seen that for a = 0, ^ is the usual fBm with Hurst parameter 
^(6 -|- 1) (up to a multiplicative constant). 

(b) Analogously to the case of fBm, in our definition of wfBm we exclude 6 = 1. It is easy to 
see that for 6 = 1 (and a > 0, according to Theorem 12. ip . p.ip is the covariance of the process 

it= I Wr'^dr, 
Jo 

where u; is a standard Bm. 

(c) For 6 = 0, ^ is a time-inhomogeneous Bm. 

In the next theorem we collect the main properties of wfBm. 

Theorem 2.4 The weighted fractional Brownian motion ^ with parameters a and b has the 
following properties: 

(1) Self- similarity: 

(^rf)t>o = (c(^+'^+')/'6)t>o for each c> 0. 

(2) Second moments of increments: for < s < t, 

E{it-is? = 2 f u'^it-ufdu, (2.3) 

J s 



A 



E{^t-Csy<c\t-sr\ (2.4) 

if a > 0, s,t < T for any T > with C = C{T), and also if a < 0, s,t > e > for any e > 0, 
with C = C{e); 

E{it-isf <C\t-s\'+'^+\ s,t>0, (2.5) 

ifa<0, 1 + a + 6 > 0; 

E{^t-Cs?>C\t-s\'+\ (2.6) 

if a > 0, s,t > e > for any e > with C = C{e) and also ifa<0, s,t <T for any T > with 
C = C{T). 

(3) Path continuity: ^ is a continuous process with the only exception of the case a < 0, b < 0, 
a + b = —1, where ^ is discontinuous at 0. 

(4) Covariance of increments: For 0<r<v<s<t, 

Q{r, V, s, t) = E{{it - is){iv - ir)) = u-[{t - uf -is- n)^]dn, (2.7) 

hence 

( >0 if b>0, 

Q{r,v,s,t) < =0 if b = 0, 

( <0 if b<0. 

(5) Long-range dependence: For r,v, s,t as in (4), 

lim T^-^Q{r, v,s + T,t + T) = — ^(t - s){v^+^ - r''+^). (2.8) 

T^oo O + 1 

(6) Asymptotic homogeneity: The finite- dimensional distributions of the process {T~"'l'^{^t+T — 
(,T))t>o converge as T ^ oo to those of fBm with Hurst paramter (1 + b)/2, multiplied by 
(2/(l + 6))V2. 

(7) Short and long-time asymptotics: 

lime-'-'E{^t+,-^tf = r^t", (2.9) 

e\0 0+1 

lim T-(i+"+'')^(et+T - 6)^ = 2 C u^il-ufdu. (2.10) 

T->oo Jq 

Hence ^ has asymptotically stationary increments for long time intervals, but not for short time 
intervals. 

(8) ^ is not a semimartingale ifb^O. 

(9) is not Markov ifb^O. 

2.2 Proofs 

Proof of Theorem 12.11 Firstly, we prove positive definiteness of Q in the case 

a>-l, -1<5<0, a + 6+l>0, (2.11) 
(see ([22])). We have from (l2T^ 

Q{s,t) = Qi{s,t)+Q2{s,t), 

where 

/•sAt t-sAt 

Qi(s,t)= / u^iis At) -u)''du, Q2{s,t) = / u^ds V t) - ufdu. 
Jo Jo 



It suffices to show that both Qi and Q2 are positive-definite. 
Qi can be written as 



Jo 



Q^{s, t) = {sA ty+^+' I u^l - ufdu, 

so it is positive definite by (|2.1ip . 

Next, since 6 < we can write Q2 as 

l-OO 

Q2is,t)= / u^[is - uf A {t - uf]l[o,s]iu)l[o,t]iu)du, 
Jo 

hence positive definiteness of Q2 follows easily. 
Now assume that 

a > -1, 0<6<1, b<l + a. (2.12) 
From ()2.ip . Q can be transformed into 

Q{s,t)=b[ I {uArf\u-r\^~^drdu, (2.13) 
Jo Jo 

hence it is clear that Q is positive-definite for 6 = 1 (note that (I2.12[) implies a > in this case). 
Assume 6 < 1; then, from (j2.13p . 

Q(s, t)=b r [\u A r)"(n V r)^~^ (l - dudr. 
Jo Jo \ u\/rj 



dudr. 



Using the power series expansion of (1 — •)'' ^ we obtain 

^/ N v^,r(l-6 + n) /■'^ f\ ^f, i(uAr)" 

Q(s,t)=y 6-^- — / / / (u Arfiuy rf-^) 

^ ' t^o r(l-6)n! Jo Jo^ ' ^ ' {uVrr 

Each summand is positive definite, since 

(u A r)"(n V rf'^u A r)"(u V r)"" = u''r''{u^-^"' A r^-i~")(n A r)"(M^" A r"") 

(we have used that 6 — 1 — a < 0, by (j2.12p ). Hence Q is positive-definite in this case. 

We will show that for the remaining values of parameters Q is not positive-definite (recall 
that Q is infinite if either a or 6 is < —1). More precisely, we will prove that there exists t > 
such that the covariance inequality 

Q(l,i) < (Q(l,l)Q(t,t))i/2 (2.14) 

does not hold. Indeed, for — 1 < 6 < 0, a > —1, a + 6 + 1 < and t \ 0, the left-hand side of 
()2.14p is of order t^+"-+'' while the right-hand side is of order t(i+a+fe)/2_ 

For a > —1, b > a + 1 and t y 00 the left-hand side of (I2.14p is of order and the right-hand 
side is of order t^^+a+b)/2^ 

It remains to consider the case 1 < 6 < a + 1. We show that (I2.14p does not hold for t = 1 + e, 
e \ 0. Using convexity of the function for 6 > 1, we have 



Q{l,l + e) > 2^ n"(^l + |-n^ 



b 

du 



£ 



a+b+l r-(l+e/2)-^ 



2(1 + -) / u%l-ufdu 



> 2( l + e+ ^ 



(a+b+l)/2 



1 



u^il-ufdu 







eV^ 1 



6+1 



This implies 

Q(l, 1 + e) - (Q(l, 1)Q(1 + e, 1 + e)f''^ 



Jo 



> 2 



2\{a+b+l)/2 N .1 

for some positive constants A, B. In the last estimate we have used the fact that a+6+1 > 26 > 2. 

The right-hand side of (j2.15p is strictly positive for e sufficiently small, since 6 > 1, so (|2.14|) 
does not hold for such e. □ 

Proof of Theorem 2.4 

(1) Self similarity follows immediately from (j2.ip . 

(2) Formula (|2.3p is a direct consequence of (|2.ip , and (|2.4p follows from ()2.3p . 
To prove (j2.5p . first observe that if a < 0, 6 > 0, (j2.3p implies 



a + 1 a + 1 

Next assume that a<0, 6<0, l-|-a-|-6>0. Fix any 9 such that |a|/(l + 6) < ^ < 1, and put 
p = 0/\a\ and q = 9/{9 — \a\). For s < t, the Holder inequality applied to (|2.3p yields 

< c(^-s)l+"+^ 

since < + 1 < 1. 

The inequality (|2.6p follows from (j2.3p . 
(3) Path continuity of ^ is a consequence of (12. 4j) and (|2.5p . If a < 0, 6 < 0, a + 6 = —1, then 



£;^2 ^ 2 /" n''(l - n)^(is for t > 0. 



(4) Formula ([IZD follows from (|2T^ . 

(5) The limit (12. 8p can be easily obtained from ([2 

(6) For < s < t, by ([l^D and 1^1^ we have 



i^i?((6+T-CT)fe+T-eT)) = fj^ u%{t + T-uf + {s + T-uf)du 



(t^+i + s^'+i _ (t _ as T^oo, 



6+1 



hence the assertion follows because ^ is Gaussian. 

(7) The proofs of ()2.9p and ()2.10p are very easy, so we omit them. 

(8) The non-semimartingale property follows from (j2.4p . (j2.6p and Corollary 2.1 in [jj. 

(9) It is not hard to see that for 6 7^ the covariance (j2.ip does not have the triangular property, 
so ^ is not Markovian (see [13], Proposition 13.7). □ 
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3 Negative sub-fractional Brownian motion 
3.1 Results 

Let us recall the function K defined by (|1.4p : 

K{s, t) = {2- h) + t'^ - + tf + \s- . (3.1) 

We know that for < /i < 2 this is the covariance of a sub-fractional Brownian motion [2j . Now 
we consider the case h > 2. 

Theorem 3.1 For 2 < h < 4 there exists a centered Gaussian process C, = {Ct)t>o with covari- 
ance 1^3. This process has a representation 

Ct = {^h{h -l){h- 2)) ' ^,ds, (3.2) 

where '& is a centered continuous Gaussian process with covariance 

Ko{s,t) = {s + tf-'^ -\s-t\^-^. (3.3) 



Definition 3.2 The process C, described in Theorem 3.1 for 2 < /i < 4 will be called negative 
sub-fractional Brownian motion (nsfBm) with parameter h. 

Remark 3.3 (a) It is not difficult to see that (|3.1|) is not positive-definite if /i > 4. Indeed, for 
t oo, K{1, t) is of order t'"''^, while K{t, i)^/^ is of order t''^'^. 

(b) For h = 4, Q = ct^^, where 7 is a standard Gaussian random variable and c is a constant. 

(c) It is easy to see by a simple covariance calculation that the process •& appearing in Theorem 
3.1 is equal in distribution to the process (^t — ^_t)t>0; where ^ is fBm with Hurst parameter 
{h — 2)/2 (defined for all —00 < t < 00). The process |(^t — i-t)t>o is called the odd part of 
fBm in [8]. 

(d) For 2 < /i < 4 the covariance of nsfBm can be written as 

[■S/\t fS ft 



K{s, t) = h{h - l){h -2Y / {u + u' - 2rY'^du'dudr. (3.4) 

Jq Jr Jr 

In the next proposition we gather basic properties of nsfBm. 

Proposition 3.4 The negative sub -fractional Brownian motion (, with parameter h has the 
following properties: 

(1) C has continuously differentiable paths, hence it is a semimartingale. 

(2) Self- similarity: 

{Cct)t>o = {c^^'^Ct)t>o for each c> 0. 

(3) Second moments of increments: for < s < t, 

EiCt - Csf < C\t - s\^ 
if s,t < T , for any T > with C = C{T); 

E{Qt-Csf>C\t-s\\ 
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"if s,t > e, for any e > with C = C{e); 

E{Ct - Csf > C\t - s\\ 

for all s,t >0. 

(4) Positive correlation of increments: For 0<r<v,0<s<t, 

mcv - Cr){Ct - Cs)) > 0. 

(5) Long-range dependence: For 0<r<v<s<t, 

lim T'-'^EiiC. - C.)(G+T - 6+t)) = _ s)iv' - r\ 

1 — >oo Zj 

(6) C is not Markov. 

We omit the proof of this proposition, as it is easily obtained from (j3.1|) and (|3.2|) . Integral 
representations of sfBm are given in [2j. Analogous results can be obtained for nsfBm from (j3.2p 
and the known integral representations of fBm (see Remark 3.3 (c)). 

Let us recall the function R defined by (|1.7p : 

R{s, t) = - (^s"^ log s + log t - ^ [(s + tf \og{s + t) + {s- tf log |s - t\]^ , s,t>0; (3.5) 

we already know that it is positive-definite. 

The proximity of the process ij = {r]t)t>o with covariance (|3.5p to nsfBm is seen if we look 
at formula (13. 4|) . since it easy to verify that ()3.5p can be written as 



R{s,t) = 2 / / {u + u -2r)-'^dududr. (3.6) 

Jo Jr Jr 

It is also worthwhile to recall that nsfBm is related to a particle system for d < a,rj cor- 
responds to the borderline case d = a, and the case a < d < 2a leads to sfBm (see the 
Introduction) . 

Some properties of rj, such as self-similarity and long-range dependence, are analogous to 
those of nsfBm and are easy to obtain. Path continuity follows for example from The 
long-range dependence is given by 

lim TE{{ri^ - T]r){r]t+T - Vs+t)) = {t - s){v'^ - r^), 

1 — ^oo 

0<r<v<s<t. The most interesting property of r] is the fact that, unlike nsfBm, r] it 
not a semimartingale, so, in this sense rj is closer to sfBm. We state this property as the next 
proposition. 

Proposition 3.5 The process r] is not a semimartingale. 
3.2 Proofs 

Proof of Theorem 3.1 Existence of the process "d follows from Remark 3.3 (c). 
From (13.11) we have 



K{s,t) = ^h{h-l){h-2) I I Ko{u,v)dudv, 



t ps 

Jo 



Q 



where Kq is given by (j3.3|) . Hence the theorem is proved. □ 

It is easy to show that Kq is not positive-definite for /i > 4, so the bfBm with H = ^ cannot 
be extended to k > 2. 

Proof of Proposition 3.5 We will show that r] has infinite variation on [0, 1], and its quadratic 
variation is zero. 

From (j3.6p for s < t, we obtain by elementary calculus 

E{rit - r]sf = I I [log(2s + u + u')- log(n - u')] du'du. (3.7) 
Jo Jo 

We puts = k/n,t = {k+l)/n,e'^ < k < n. Using log (2A;/n) < log{2k /n+u+u') < log((2A;+2)/n) 
for < ti, < 1/n, it is not hard to see that ()3.7p implies 

^ log A: < E{r](^k+i)/n - Vk/nf < ^ log k. (3.8) 

(The assumption A; > is used in the upper estimate, where we need to have that (log (2 A; + 
2) + 2)/2 < log A;. The lower estimate holds for any k > 1.) 
From the upper estimate in (j3.8p we infer that 

n-1 



Jl™^ ^ E{r](k+l)/n - Vk/nf = 0, 



fc=0 

hence the quadratic variation of r] cannot be different from 0. 
By the Gaussian property of r/, (|3.8p implies that 

(log A;) 

E\ri(k+l)/n - %/nl > \^ ■ (3.9) 

To show that the variation of ry is infinite, we repeat the argument of the proof of Lemma 
2.1 in [Ij. By the — 1 law for Gaussian processes, the variation on [0, 1] is either finite a.s. 
or infinite a.s. If it were finite a.s. (we denote it by |r/|i), then by the Fernique theorem [9] we 
would have E\r]i\i < oo. On the other hand (j3.9p implies that 



n-1 

,\1/9 

oo as n — > oo. 



^ n— 1 

E\r,U>-Y,{logk)'/' 

m ^ — ^ 



n 

k=l 



Therefore I77I1 = 00 a.s. □ 
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sity of Mexico (UNAM), where this paper was written. 
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